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Abstract 



One of the main applications in plasma physics concerns the energy production through 
thermo- nuclear fusion. The controlled fusion is achieved by magnetic confinement i.e., the 
plasma is confined into a toroidal domain (tokamak) under the action of huge magnetic 
fields. Several models exist for describing the evolution of strongly magnetized plasmas, 
most of them by neglecting the collisions between particles. The subject matter of this 
p ^ ' ' paper is to investigate the effect of large magnetic fields with respect to a collision mech- 

anism. We consider here linear collision Boltzmann operators and derive, by averaging 
with respect to the fast cyclotronic motion due to strong magnetic forces, their effective 
i-^ , collision kernels. 
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■ 1 Introduction 

(y-^ ' The goal of this paper is to emphasize the role of the collision mechanism in the framework 

of magnetic confinement. We investigate both first and second order asymptotic approxima- 
tions, when the magnetic field becomes large. We present completely explicit models and 
clearly indicate the form of the effective collision operator, the expression of the collision 
drift and check that the standard physical properties hold true. 

Motivated by the magnetic fusion, many research programs in plasma physics concern the 
dynamics of charged particles under the action of strong magnetic fields |B £ | — > +oo as e \ 0. 
Usually we appeal to kinetic models meaning that we replace the particles by their density 
(distribution) function f £ = f £ (t,x,v) over the phase space, depending on time t G M+, 
position x G 1R 3 and velocity v G M 3 . The quantity f £ (t,x,v) dxdv represents the particle 
number at time t, inside the infinitesimal volume dxdv around (x, v). The fluctuations of the 
density f £ are due to the transport of the charged particles under the electro-magnetic field 
and also to the collisions between particles 

d t f + v-V x f + ^-(E(x)+vAB £ (x))-V v f £ = Q(f £ ), (t, x, v) G R + x M 3 x M 3 . (1) 
m 
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We add the initial condition 



f(0,x,v) = f in (x,v), (x,v)€R 3 xR 3 . (2) 

Here m is the particle mass and q is the particle charge, (E, B e ) stands for the electro- magnetic 
field and Q is the collision operator. We neglect the self-consistent electro-magnetic field; we 
assume that the magnetic field is stationary, divergence free and that the electric field derives 
from a given electric potential E(x) = —V x (p(x). We focus on the linear Boltzmann collision 
operator 

Q(f)(v) = - I s(v, v'){M(v)f(v') - M(v')f(v)} dv' (3) 



r _ 

where M is the Maxwellian equilibrium 

2 



M(v) = — —j- exp 



m\v\ 



(27r0/m) 3 / 2 y 26 

with temperature 9 > 0, r is the relaxation time and s is the scattering cross section. We 
assume that 

s(v,v') = s(v',v), v,v' e R 3 (4) 

and 

so < s(v,v') < S , v,v' e M 3 (5) 

for some so,So > 0. We analyze the asymptotic behaviour of ([1]) when the magnetic field 
becomes large 

TS(x) 

B e (x) = — B(x) = B(x)b(x), div x (Bb) = 0, < e « 1 (6) 

for some scalar positive function B(x) and some field of unitary vectors b(x). 

These studies are motivated by the numerical simulations of strongly magnetized plasmas. 
Indeed, a direct resolution of ([T]) using explicit numerical schemes requires a time discretiza- 
tion of order C(l/|B e |) = 0(e), and therefore a huge computation time, when e becomes 
very small. The main idea will be to perform the asymptotic analysis of ([1]) at the theoretical 
level i.e., to replace ([T]) by first/second order approximation models whose coefficients are 
stable when e \ 0. Eventually we solve numerically these asymptotic models. Notice that 
in this case the largest time step satisfying the stability condition does not depend on e, and 
therefore the numerical resolution remains efficient even for very large magnetic fields. 

For simplicity we neglect the magnetic curvature i.e., d x b = 0. Assuming that b = e% 
and taking into account the divergence constraint in ([6|) lead to B e (x) = (0, 0, B £ (x\, Xy)) = 
(0, 0, B{x\, X2)/e). Therefore the equation (P) should be understood in two dimensional 
setting 

dtf + v ■ V x f + ± (e(x) + —B(x)) ■ V v f = Q(f), t > 0, (x, v) e R 2 x M 2 (7) 

where v = (t>2, — i>i) and Q is the two dimensional linear Boltzmann operator, associated to 
the Maxwellian M(v) = m/(2vr(9) exp(-^f),v G R 2 . 

The dynamics of the particles is dominated by the transport operator where 

r=^M^-v„. (8) 

m 
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When e is small enough it is natural to introduce expansions like f £ = f + sf 1 + £ 2 / 2 + ••• 
Plugging the previous Ansatz in ([7]) gives at the lowest order the divergence constraint Tf = 
and to the next order the evolution equation 

dtf + v- V x f + ±E(x) ■ V v f + Tf 1 = Q(f). (9) 
m 

We expect that the Boltzmann equation ([9]) supplemented by the constraint Tf = (and 
also by an appropriate initial condition) provides a well posed problem for the leading order 
term /. 

The approach is standard as in the Hilbert expansion methods of matching orders of e: 
the first order term f 1 enters |9| as a Lagrange multiplier in such a way that at any time 
t > the constraint Tf(t) = should be satisfied. Finding a closure for the dominant term / 
reduces to eliminating the Lagrange multiplier provided that the constraint on / holds true. 
Notice that the range and kernel of the operator T are orthogonal (in L 2 ) and therefore the 
closure for / follows naturally by taking the orthogonal projection on kerT of the Boltzmann 
equation at any fixed time t > 0. One of the key point is to observe that taking the 
orthogonal projection on the kernel of T reduces to averaging along the flow associated to 
T cf. [3], [I], [5]. Observe that the motion associated to the vector field (0, 0, uj c (x) ± v) is 
periodic 

X(s;x,v)=x, V(s; x,v) = R(—uj c (x)s)v 

where oj c (x) = qB(x)/m is the rescaled cyclotronic frequency and for any a£l the notation 
R(a) stands for the rotation of angle a 



R(a) 



cos a — sin a 
sin a cos a 



We introduce the rescaled cyclotronic period T c (x) = 2tt/lu c (x). In this case, for any L 2 
function u, the average simply writes 

1 fT c (x) 

(u)(x,v) = / u(X(s; x, v), V(s; x, v)) ds 



T c (x) jo 



= — / u(x, R(— lo c (x)s)v) ds. 
2?r Jo 

Applying the average operator to the Boltzmann equation (|9|) yields the following closure for 
the dominant term / 

(d t f) + (v ■ V x f + S- E {x) ■ V v f) = (Q(f)) (10) 

supplemented by the constraint Tfit) = at any time t > 0. Notice that {ip\ = xi,ip2 = 
x 2, ip3 = Ml form a complete family of prime integrals for T and thus the constraint Tf = 
reads f(t,x,v) = g(t,x,r = \v\) for some function g = g(t,x,r). Straightforward computa- 
tions show that 

m) - a, (/, = o tf , (, • v,/, = („ ■ v, s(t , „ H) > — o, (e- v„/> = (e . ^a„) = o. 

Moreover, we will see that the linear collision operator commutes with the average operator 
i.e., {Q(f)) = Q({f)) = Q(f) and finally the dominant term satisfies the space homogeneous 
Boltzmann equation |12| . [15j . |25j 

<hf Q(f). 
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In particular there is no current j R2 vf(t,x,v) dv = f R2 vg(t, x, \v\) dv = along the or- 
thogonal directions to the magnetic lines, as expected by the magnetic confinement context. 
But orthogonal drifts may appear at the next order. It is well known that the particles 
move slowly along the orthogonal directions and we distinguish the electric cross field drift 
v £ h = the magnetic gradient drift ^q D = — i^j— tjf^-,^ = the magnetic curvature 
drift (which is neglected here since d x b = 0) cf. [27]. We intend to investigate the effect of 
strong magnetic fields through the collision mechanism. We will also prove that the collisions 
account for a drift motion which, to the best of our knowledge, has not been reported yet. 
In the case of a constant scattering cross section s we find the formula 

e s ±v 

V QD = 7 

(see (|69p for non constant scattering cross section). Naturally, in order to observe the previous 
drifts (of order of e) we need to take into account the first order correction term in the 
expansion f 6 = f + ef 1 + e 2 f 2 + ... eventually to appeal to the next order evolution equation 

dtf 1 + v ■ V,/ 1 + ?-E ■ V v f l + Tf = Q(/ 1 ). (11) 
m 

As before we are looking for a closure with respect to / + ef , by getting rid of the second 
order correction f 2 . This will be achieved using again the averaging technique and combining 
properly the time evolutions of /, f 1 in @, (|11|) . Doing that we will get another Boltzmann 
like equation, whose collision operator Q + eQ 1 is a first order perturbation of the standard 
linear Boltzmann operator Q. The interesting point is that the new collision operator is non 
local in space anymore. For example, when the scattering cross section is constant we obtain 

eQ\f)(x,v) = -M(v) im ^ p f (x) - — £ ■ L V xPf + diy 

where pf(x) = L 2 /(a?,?/) dv',jf(x) = f R2 v'f(x, v') dv' . Surely the construction of the 
new collision operator should take into account the classical physical properties, as mass 
conservation and entropy inequality . We will see that eQ 1 (and therefore Q + eQ 1 as well) 
satisfies global mass conservation and entropy inequality. 

We point out that a linearized and gyroaveraged collision operator has been written 
in [29], but the implementation of this operator seems very hard. We refer to [10J for a 
general guiding-center bilinear Fokker-Planck collision operator. Another difficulty lies in the 
relaxation of the distribution function towards a local Maxwellian equilibrium. Most of the 
available model operators, in particular those which are linearized near a Maxwellian, are 
missing this property. Very recently a set of model collision operators has been obtained in 
[T7j . based on entropy variational principles. 

The analysis of the Vlasov or Vlasov-Poisson equations with large external magnetic 
field have been carried out in [IB], |2Uj . The numerical approximation of the gyrokinetic 
models has been performed in [21] using semi-Lagrangian schemes, see also [15], |19| . For 
recent results on the Boltzmann equation we refer to [7] , [8J . High electric field limits of the 
Boltzmann-Poisson system have been addressed in [26], [13], |14| . pp. 

The nonlinear gyrokinetic theory of the Vlasov-Maxwell equations can be carried out by 
appealing to Lagrangian and Hamiltonian methods [11] , [23] , [23] . It is also possible to follow 
the general method of multiple time scale or averaging perturbation developped in [2j. For a 
unified treatment of the main physical ideas and theoretical methods that have emerged on 
magnetic plasma confinement we refer to [22] . 

The outline of the paper is the following. In Section [5] we introduce our models and 
present the main results. Section [3] is devoted to the definition and main properties of the 
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average operator. In particular we investigate its commutation properties with respect to first 
order differential operators and the collision operator. The first and second order asymptotic 
models for the linear Boltzmann equation with large magnetic field are detailed in Section [H 



2 Presentation of the models and main results 

As observed formally before, the first order approximation of ([7J) leads to the space 
homogeneous Boltzmann equation. We obtain rigorous weak and strong convergence results, 
at least for electric potential satisfying 

E{x) = -V x <j){x), J exp(-q^^\ dx < +00. (12) 

The equilibrium F(x,v) = M{v) exp (—q<p(x)/Q) belongs to L 2 F = L 2 (-^rfy) and by stan- 
dard computations we obtain the uniform bound 

<>■<'■*■•»'**«< / / d»d,, ,>0 



F(x,v) ' ' Jr2J R 2 F(x,v 

which allows us to get weak compactness in L°°(1R + ; L 2 F ). For scattering cross section of the 
form s{v,v') = a(\v — v'\) one has 

Theorem 2.1 Assume that the scattering cross section has the form s(v,v') = o~(\v — v'\) G 
[so, Sq], v, v' G R 2 and that the electrostatic potential verifies (|12p . For any e > we denote 
by f e the weak solution of (|7|) , ([2|) associated to a initial condition satisfying f m G Lp. 
Therefore (/ e ) e >o converges weakly * in L°°(IR + ;L^) as e \ towards the solution of the 
problem 

d t f = Q(f), (t, x, v) G R\ x R 2 x R 2 , /(0, x, v) = </ in > (x, v), (x, v) G R 2 x R 2 . (13) 

Strong convergence results hold true as well. The main idea is to introduce a smooth first 
order corrector, under appropriate smoothness hypotheses on the initial condition. 

Theorem 2.2 Assume that the scattering cross section satisfies s(v,v') = a(\v — v'\) G 
[so, So], v, v' G R 2 for some function a G W 2, °° , that the electro-magnetic field (E(x) = 
-V x (/),B(x)) belongs to W 1 ' 0O (R 2 ) such that 

[ exp ( ) dx < +00, inf B(x) > 

and that the initial condition verifies 

rr = o, (i + \v\ 2 )\vi v r\eL 2 F . 

Then for any T > there is a constant Ct such that 

\\f(t)-f(t)\\ L 2 <C T e, te[0,T], e>0. 

F 

The leading order term in the expansion of f £ satisfies, at any fixed space point, the homo- 
geneous Boltzmann equation. There is no phase space transport, the collision operator Q is 
local in space and there is no current, since f(t,x, •) depends only on \v\. A much more in- 
teresting model will be the second order approximation, since we expect to retrieve the usual 
drifts along the orthogonal directions with respect to the magnetic lines, but also a collision 
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drift. Certainly, the second order approximation will appear as a first order perturbation of 
the space homogeneous Boltzmann equation, involving a new transport operator but also a 
new collision operator Q + eQ 1 , which is non local in space anymore. The full expression 
of the operator Q 1 is rather complicated, see Corollary 14.11 but it simplifies a lot when the 
scattering cross section is constant 



Q\f) = -M{v) 

T 



oj c (x) J u c (x) \Pf(x 



where Pf = L2/ dv', pp = dv', jf = J^v 1 f dv'. For the sake of the presentation 

we formulate the second order convergence result in the context of constant scattering cross 
section but we point out that similar results hold true in the general case. 

Theorem 2.3 Assume that the electro-magnetic field {E(x) = —V x (j){x),B £ (x) = B{x)/e) 
is smooth such that 



f 

Jr 



exp ( -^M ) dx < +00, inf B(x) > 0. 



We denote by f £ the weak solution of ([7]) with the initial condition f £ (0) = f[ n and by f e the 
weak solution of 

„f £ ^f ± E \v\ 2 ^VJ^\ ~ ( ± E V X B £ \ ~ 



\B £ 2oo e c B £ J J \2B e B 
with the initial condition 

f(o) = r + -4 • v,/ in - -J • v,/ in , rr = (15) 

u>l B £ 

where uj £ c = 2— = ^ is the cyclotronic frequency. If the family of the initial conditions 
(/in)e>o satisfies 

fL = r + ^ • v,r - ^ • v.r + o(e 2 ) 

then we have f £ = f e + 0(e 2 ). 

Remark 2.1 Notice that the coefficients in (|14p still depend on e but they are stable when 
e \ 0. In particular the numerical approximation can be achieved uniformly with respect to 
e > 0, using a time step not depending on e. 



3 Averaging along a characteristic flow 

This section is devoted to the average technique. It consists in averaging along the 
characteristic flow associated to a divergence free vector field. When applied to the gyrokinetic 
theory, this method is known as the gyroaverage technique. It allows us to derive limit models 
by neglecting the fluctuations associated to fast motions, as the cyclotron rotation of particles 
around the magnetic lines. From the mathematical point of view, it reduces to averaging along 
the characteristic flow associated to the dominant advection field (here the magnetic force 
in the framework of magnetic confinement) and it is also equivalent to projecting onto the 
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subspace of functions which are left invariant by the same flow. These tools are well known 
and have been used when establishing the guiding center approximation in three dimensions 
for general magnetic shapes [5], the finite Larmor radius regime [3] and, more generally, 
when studying transport equations with disparate advection fields [IJ. We point out that 
the averaging technique applies to many other cases : models which take into account the 
mass ratio between electrons and ions subject to strong magnetic fields [9] and models with 
fast oscillating magnetic fields [6]. In the present work we generalize the averaging method 
to collisional models. Motivated by the study of the Boltzmann equation, we revisit some 
results concerning the average of first order differential operators which will combine with 
other new results in order to study the average of collision operators. 



3.1 Gyroverage operator 

We recall the definition and the properties of the average operator corresponding to the 
transport operator T given in ([5]), now in divergence form, whose definition in the L 2 {R 2 xR 2 ) 
setting is given by 

_ ,./,., . I \ . . qBix) 

Tu = div v uiAx) u{x, v) v , ojAx) = 

V / m 

for any function u in the domain 

D(T) = {u(x,v) G L 2 {R 2 x R 2 ) : div,, (u c (x) u ± v) G L 2 (R 2 x R 2 )}. 

We denote by || • || the standard norm of L 2 (R 2 x R 2 ) and by (X, V)(s; x, v) the characteristics 
associated to ui c {x) v ■ V v 

^ = 0, ^= U ; c (X(s)) ± V(s), (X,V)(Q) = (x,v). (16) 
as as 

By direct computations we deduce that 

X(s;x,v) = x, V(s; x,v) = R(—oj c (x)s)v. (17) 

A complete family of invariants for the flow (TT7|) is given by {21,22, M = ((^l) 2 + (v2) 2 ) 1 ^ 2 } 
and therefore the kernel of T is given by those functions in L 2 (R 2 x IR 2 ) depending only on 
21,22 and \v\. Notice that all trajectories are T c (x) = 2it/uj c (x) periodic and we introduce 
the average operator cf. [4J 

1 fTc(x) 

(u)(x,v) = 7f~~r^ / u(X(s;x,v),V(s;x,v)) ds 
T c {x) Jo 

Y r2ir ^ r 2ir 



u(x, R(—a)v) da = — / u(x, R(a)v) da 

for any function u G L 2 (R 2 xR 2 ). The following two results are justified in Propositions 
2.1, 2.2. The first one states that averaging reduces to orthogonal projection onto the kernel 
of T . The second one concerns the invertibility of T on the subspace of zero average functions 
and establishes a Poincare inequality. 

Proposition 3.1 The average operator is linear continuous. Moreover it coincides with the 
orthogonal projection on the kernel ofT i.e., 



<«> 



G kerT and / / (u — (u))(p dvdx = 0, V(/?Gker7". 
Jr 2 Jr 2 
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Remark 3.1 Notice that we have the orthogonal decomposition of L 2 (M 2 xl 2 ) into invariant 
functions along the characteristics (|17p and zero average functions 

u = (u) + (u — (u)), / (u — (u)) (u) dvdx = 0. 
Jr 2 Jm? 

Taking into account that T* = —T, the equality (■) = Proj ker7 - implies 

ker (•) = (kerT)- 1 = (kerT*)- 1 = Range T. 
In particular Range T C ker (•). 



The next point is very important as it relates to the solvability of Tu = v. We prove that 
the range of T is closed and derive the corresponding Poincare type inequality, provided the 
magnetic field remains away from 0. 

Proposition 3.2 We assume that inf^g^ B{x) > 0. Then T restricted to ker (•) is one 
to one map onto ker (•). Its inverse belongs to £(ker (•) , ker (•)) and we have the Poincare 
inequality 

\\u\\ < | — r||Tti||, ojq = — inf B(x) ^ 

| wo I fn x<m? 

for any u £ D(T) fl ker (•) . 
Another useful result is given by 

Proposition 3.3 Assume that (p is a C 1 (1R 2 x R 2 ) function such that ip and Tip are bounded 
on M 2 xl 2 . Then, for any function u G D(T) we have the integration by parts formula 

(pTu) + (uTip) = 0. 

Proof. We have ipu G D(T) and T(<pu) = ip Tu + u Tip. Taking the average yields 

{tp Tu) + {u Tip) = (T((pu)) = 

since Range T C ker (•). |— I 



3.2 Averaging and first order differential operators 

In the sequel we need to analyze the commutation properties of averaging and transport 
operators. For the sake of the presentation we perform our analysis in the general framework 
of characteristic flows associated to divergence free vector fields. We denote by b° a given 
smooth, divergence free vector field b° : R m — > W 71 whose characteristic flow Y = Y(s; y) is 
well defined 

/ £Y =b°(Y(s;y)), (s,y)eMxR m 

\ Y(0;y) =y, y G R m . 1 °> 

We intend to apply the following general results to the gyroaverage context and therefore, 
without loss of generality, we will assume that the characteristic flow Y = Y(s; y) is periodic 
i.e., for any y G M. m there is T c (y) > such that 

Y(T c (y); y )=y. (19) 

We may assume that T c (y) is the smallest positive period, excepted at those points y G M m 
where b°(y) = when s — > Y(s; y) is constant. The non periodic case is beyond the scope of 
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this work. Nevertheless notice that it can be handled as well, the average operator associated 
to a non periodic flow being defined through ergodic mean (see [I] for details) 



1 f T 

(u){y)=l™ 7f u(Y(s;y)) ds, yG 



For the sake of clarity, the results formulated in the general framework are stated by using 
latin upper case T,A,B,C,... whereas their counterpart results, referring to the gyroaverage 
framework, are stated by using calligraphy upper case T, A, B,C, .... 

Proposition 3.4 The map y — > T c (y) is constant along the flow Y . 

Proof. Indeed, for any y G M m such that b°(y) ^ let us consider y^ = Y(h; y), with fteR. 
Clearly b°(yh) ^ (otherwise Y(-; y) = y and b°(y) = b°(yh) = 0) and 

Y(T c (y); y h ) = Y(T c (y) + h;y)= Y(h; Y(T c (y); y)) = Y(h; y) = y h , 

which implies T c (y/j) < T c {y). Replacing now y by y„h (which satisfies 6°(y_/j) / 0) we 
deduce that T c (y) < T c (y_^) for any h € R and therefore, after changing h by — h in the last 
inequality we obtain T c (yh) < T c (y) < T^y^), saying that T c (y) is constant along the flow Y. 

□ 

Next, we denote by T the linear operator defined by Tu = div y (u(y)b° (y)) for any u in 
the domain 

D(T) = {u € L 2 (R m ) : div y (b° (y)u(y)) G L 2 (R m )} , 
and thus the kernel of T is given by constant functions along the characteristic flow ()18p 

kerT = {u G L 2 (R m ) : u(Y(s;y)) = u(y), s G R, a.e. y G R m }. (20) 

For further computations it is very useful to pick appropriate coordinates. Since we are 

interested in averaging along a characteristic flow (and therefore since the functions left 
invariant along this flow play a crucial role) we take as new coordinates a complete family of 
invariants {ipi, ...,ip m -i} for the field b°, together with a parametrization tpQ along its flow. To 
any invariant ipi, i G {1, m— 1} we associate a derivation b l ■ V y acting, with respect to the 
new coordinates, like the partial derivative . Assume that the vector field b° has a complete 
family of smooth independent prime integrals denoted ipi = ipi(y), ...,ip m -i = ip m -i(y) i.e., 

b°-V y tpi = 0, i G {1, m — 1}, rank3/(^i,...,Vm-i) = m-l, y G R m . (21) 

Moreover, assume that there is a function tpQ = ipo(y), whose gradient V y ipo is well defined 
a.e. y G R m , such that 

b° ■ V^o = I(y) + 0, a.e. yGR m (22) 

for some function / G ker T. Notice that if>o, satisfying (|22p . cannot be continuous everywhere. 
In general, tpo has a discontinuity point, say yo, on each characteristic and the jump [ipo] 
around the discontinuity is given by 

lim MY(s;yo)) - lunMY(s;y )) = (b° ■ V y ^ )(Y(s; y )) ds = T c (y )I(y ). 

(23) 

That means ipo(y) is an angular coordinate. Since T C (-),I(-) G ker(6° • V y ) it is possible to 
normalize the jump. Indeed, taking ipo = ipo/[ipo] leads to 
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Therefore we assume, without loss of generality, that the angular coordinate ■i/'o has the 
same jump on every characteristic. The following result is standard. For the sake of the 
completeness we detail the proof in Appendix [Aj 

Proposition 3.5 Assume that the vector field b° has a complete family of smooth indepen- 
dent prime integrals (|2ip . Moreover we assume that there is an angular coordinate (|22p 
with constant jump S = T C I on every characteristic. Then, the following two properties are 
satisfied 

1. The map y — > (ipo(y), ipi(y), ipm-i (y)) is a bijection between M m and [0,5) x D, 
where D := {(Mv), -^m-i(y)) ■ y G M m }. 

2. For any i G {1, ...,m — 1} there is a unique vector field b % such that 

b i -V y ^ j = 5), je{0,l,..,m-l}. 

Remark 3.2 Any function u = u(y) writes u(y) = w(tpo(y),ip\(y), ...,ip m -i(y)), with w de- 
fined on[0,S)xD. Notice that ifu is continuous then \\m. r/ *s uj(r,ipi, ip m -l) = w(0,tpi, ?/> m _i). 
Indeed, for any (ipx, ...,ip m -i) G D let us denote by yo the discontinuity point of ipo along the 
characteristic associated to the invariants (ipx, ip m -i). We have 

tu(0,^i,...,^ m _i) = u(y )= lim u(Y(s;y )) 

s/*T c (y Q ) 

= lim w(^ {Y{s;y )),^i{yo),...,ip m -i(yo)) 
= lim w(ip (y )+sI,Tp 1 ,...,iJ; m -i) 

= ]imw(r,ipi,...,ij) m -i). 

r y-s 

If u is continuous, the function w above extends by S- periodicity with respect to ipo to a 
continuous function, still denoted w, defined onM x D. 

Remark 3.3 In our case ([7]) the divergence free vector field to be considered is b° : K — > Mr, 
b°(x,v) = (0, 0, u c (x)- L v). A complete family of prime integrals for this vector field is 

V>1 = Xi, 1p 2 = X 2 , 1p3 = \v\ 

and we pick as angular coordinate ipo(x, v) = —a(v) where a(v) is the angle ofve M 2 \{(0,0)} 



i.e., 



v% = \v\ cos ct(v), V2 = \v\sina(v), a(v) G [0, 2ir). 



The angular coordinate ipo(x,v) is discontinuous on M+ x {0}. But its gradient, which is well 
defined for any v £M 2 \ (M + x {0}), can be extended by continuity on M 2 \ {(0, 0)} 



v I v 

V v a = --—y, V XjV ip = 0, 0, r-ry 
\v\ V \v\ 



and we have I = ui c (x) v ■ V v ipo = w c (x) G ker T. On each characteristic ipo has a jump of 
S = T c (x)ui c (x) = 2tt. 

In addition, the vector fields (& l )ie{i ) 2,3} associated to the invariants ip\ = x\,ip2 = 
%2,ip3 = \v\ and the angular coordinate ipo = —a(v) are given by 

b x = (1,0,0,0), b 2 = (0,1,0,0), b 3 =(o,o,^-Y 
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The vector fields (6 J ) ie { lj ... jm _n will play a major role in our analysis cf. [5]. We shall 
establish several general results which mainly concern the commutation properties between 
the average and divergence operators, the inversion of T = b° • V y on zero average functions 
and the construction of some auxiliary vector fields we will need for our further analysis. 
In particular we will show that after averaging, a transport operator reduces to another 
transport operator and we explicitely compute its advection field components. The following 
eight propositions explain in detail how to achieve these goals. Observe that the derivations 
along the vector fields (&*)ie{i coincide with the partial derivatives with respect to the 

invariants (^i)ie{i,...,m-i}; see Appendix lAl for details. 

Proposition 3.6 Consider a function u = u(y) and let us denote by w = vj(i/jq, ip m -i) : 
lxD->M the S periodic function with respect to ipo such that u(y) = w(ifjo(y), ...,ip m -i(y)). 
If u G C^R" 1 ), then w G C 1 (R x D) and satisfies the identities 

b° ■ V y u(y) = I(y) d^ w(ip (y), ...,if> m -x(y)), V ■ V y u{y) = d$.w{i> Q (y), i> m -l {v)) 

for any i € {1, m — 1}. 

In the sequel, for any two smooth vector fields b(y) and c(y), the notation [b, c] stands for 
the vector field (b • V y )c — (c • V y )b, which is the vector field of the transport operator 
b • V y (c • V„) - c • V y (b • V,,) 

b ■ Vy( C ■ Vy) - C • Vy(b ■ Vy) = [b, C] • Vy. 

The commutators between the fields (£>*)ie{o,...,m-l} come easily by Schwarz's theorem on 
mixed partial derivatives [28J. 

Proposition 3.7 With the previous notations we have 

[ b \ b o ] = l^lL b O ) i g{ ,...,m-l} 
and [&,b>] = 0,i,j £ {1, ...,m - 1}, with I = b° ■ V y ip . 

Proof. The first statement is obvious when i = 0. For any smooth function u = u(y) and 
i £ {1, ...,m — 1} we have, after the change of coordinate u(y) = w(ip(y)),ip = (ipo, ip m —i) 

[b\ b ] ■ VyU = (b l ■ Vy)(b ■ VyU) - (b ■ Vy){U ■ VyU) 

_ (6 -.v s ,(^o,)- (6 0. Vs) (^ o , 

= v — b° ■ V v u. 

Similarly, for any i,j G {l,...,m — 1} we have 

[b\ b>] ■ V y u = (b l ■ V y ){V> ■ Vyu) - Qp ■ V J/ )(6 i • V y u) 
d ( dw \ \ , ( d ( dw \ \ 



dipi \ difjj J J \ dijjj \ dij). 

□ 

Corollary 3.1 The operators (b l ■ Vj)^ ^.u leave invariant ker6° • V^. 
Proof. By Proposition 13.71 we have for any u G ker(6° • Vy) Pi D(b l ■ Vy) 

(V ■ Vy)(6 • VyU) ~ (b • Vy)(W ' V„«) = (6° ■ V yU) 

and therefore (6° • V y )(6* • V y u) = 0. □ 
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We establish now some properties regarding the divergences of the fields 

Proposition 3.8 The matrix (b l ■ V J/ (div ?/ 6 J ))(j J ) e {o,...,m-i} * s symmetric. In particular, for 
any i G {1, ...,m — 1} the divergence of b' L is constant along the flow ofb . 

Proof. By Proposition 13,71 we know that 

[b\V} = \ ij b°, b°-V y X ij = 0, i,j G {0, m — 1} 

since I and b l • S7 y I are constant along the flow of b° cf. Corollary 13.11 Therefore we have 

= div y (A«6°) = div„[&*, &»'] = {b l ■ V y )(div y V) - (b> ■ VyXdiVyb*). 

In particular, taking j = and i G {l,...,m — 1} one gets 6° • V J/ (div 2/ 6 i ) = 0, since b is 
divergence free. ^ 

We investigate now the commutation properties between the operators (6 l -Vy)j 6 / 0) ... )m _i} and 
the average operator. Since the left hand side of the Boltzmann equation (|7|) can be written 
into conservation form, it is worth analyzing how to average the divergence of a vector field. 

Proposition 3.9 With the previous notations, we have for any smooth vector field £(y) = 

(£i(y),-,£ m (y)) 

(div tf = div, |E (C • V v fc) 6 l | = div„ | + E ^ ■ V ^»> hi 

Proof. Notice that div^ | ftb.y"^ &° | = and therefore it is enough to prove only the first 

equality. Obviously, the left hand side of the first equality is constant along the flow of b°. 
The same happens for the right hand side since 

{ra— 1 ^ m— 1 m— 1 

E • v ^> 6 7 = E ^ • div ^ + E fei • v v • 
i=i j i=i i=i 

and we know that div y b l G ker(6°-V y ) and that (b l - V^g/i , jm _n leave invariant ker(6°- Vy). 
Therefore, in order to prove the first equality, it is sufficient to test against smooth functions 
ip G ker(6° • V^). By the definit ion of the fields (6*)ig{i ) „. )m _i} we have 

£ = ^^° + eVv^ 

i=\ 

and thus for any smooth function if G ker(6° • V y ) we obtain 

m— 1 

£-V^= E(^- v ^)(^-V^). 

i=l 
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After integration by parts one gets (by taking into account that ip, b l ■ V y ip G ker(6° • V^)) 
f(y) ( div y0 dy = / ip(y)div y £ dy 



- / £-V y </?dy 

m—1 . 

-E / (e-V^)(^-V^)dy 

i=l ■ /Rm 

m—1 „ 

-E / (e-V^^-V^dy 

i=l 

/ </?(y)div y J V (£ ■ V y Vi) 1 dy 
m I i=i J 



saying that 



(m-l \ 



□ 



Corollary 3.2 Assume that the fields (&*)iefl,...,m— 1} are smooth and have bounded diver- 
gence. Then the operators (b l ■ Vy)ie.ro m — n are commuting with the average operator. 

Proof. Obviously 6° • V y and (•) are commuting, since for any function u G D(6° • V^) we 
have 

(b° ■ V y u) =b° • V„ («) = 0. 

For any i G {1, ...,m — 1} and function u G D(6* • Vy) we obtain thanks to Proposition 13.91 
applied to the field ub l 

(div„(u&*)) = div„((u) 6*) 
and our conclusion follows easily, since we know that d\v y b l G ker(6° • V y ). | — | 

Remark 3.4 7n the particular case of the vector fields b° = (Q,0,u c (x) v), 0,0,0), 
b 2 = (0,1,0,0), b 3 = (0,0,v/\v\) we obtain for any smooth vector field £ = (£r(x, v), £ v (x, v)) 
the formulae, thanks to the identity div v {v/\v\ 2 } = 

v 

(div x>v £) = div x (£ x ) + • V„ (£„ • v) , (div^) = (div x>v (£ x , 0)) = div x . 

The previous result guarantees that a transport operator reduces, after average, to another 
transport operator and allows us to keep the model into conservative form cf. [5]. For the 
sake of transparency we detail this computation in the following proposition. 

Proposition 3.10 Consider a = a(y) a smooth divergence free vector field in R m . Then 
there is a divergence free vector field A = A(y) such that for any smooth function u G 
ker(6° • V„) U {Vo} 

(div y (ua) )= div y (u A). (24) 



13 



Proof. We apply Proposition 13.91 to the field y — > u{y)a{y) 



<div tf (ua)> = div, j (UQ " J^ 0) b° + ( ua ■ V ^> &' j • 
If u G ker(6° • V y ) then 

(div„(ua)) = div y ^ J a ' V ^°\ o + u (a • V^i) b* \ . (25) 
In particular for u = 1 one gets 

r (a-V^o), "'~ J ^ 



m— 1 



= (div„(a)) = div, y a ' y W b° + E {« • V,^) 6* }- . (26) 
Hence, we define the field A by 



i=l 



A= (a-V^ ) b o + g (Q , v ^. ) y (27) 

1=1 



so that ^4 is divergence free and the equation (|25j) becomes 

(dw y (ua)) = div y (uA), u G ker(6° • V y ). 

It remains to check that the same equality holds true for the angular coordinate tpQ. Indeed, 
we have 

(div y (V>oa)) = (a ■ V y 4> ) = A ■ V^Vo = div y (V>o-4)- 

Actually the vector field A constructed above is the only vector field satisfying (|24p . since 
necessarily satisfies A ■ V y ipi = (a ■ Vytpi) for any i G {0, m — 1}. | — | 

Remark 3.5 If {olq, a m _i) are the coordinates of a with respect to the basis {6°, b" 1 ^ 1 }, 
then ((«o) > ( a m— i)) are ^ e coordinates of A with respect to the same basis 



TO— 1 TO— 1 

a = 

i=0 i=0 



Remark 3.6 /n i/ie particular case of the vector fields b° = (0,0, lo c ( x ) v), fe 1 (l, 0,0,0), 
6 2 = (0,1,0,0), 6 3 = (0,0, f/H), a= (v,q/mE) we obtain the coefficients 



qE v qE v 

TT ' TT? ' a i =v ii a 2 = v 2 , a 3 = — 

muj c {x) \v\* m \v\ 



It is easily seen that (oj) = 0, i G {0, ...,3} and f/ms .4 = Yli=o ( Q «) ^ = 0- 



Notice that for any i G {0, ...,m — 1} the function aj — has zero average, and thus, by 
Proposition 13. 2\ there is a unique satisfying 6° - V y f3i = on — (en), with (/%) = 0. Therefore, 
we denote by B the vector field 

TO— 1 

5=£ ftfe* (28) 

8=0 

whose characterization makes the object of the next proposition. 
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Proposition 3.11 Consider a = a(y) a smooth, divergence free vector field in K m . Let A 
and B be the vector fields defined in (|27p and (|28|) respectively. 

i) Then for any smooth function u G ker(6° • V«) U {ipo} we have 

div y {(B ■ V y u)b } = div y {u(a - A)} (29) 

and diVy(IB) = 0. 

ii) Let B = B + X b° where b° ■ V y A = ^x± } (\ ) = 0. Then for any smooth function 
u G ker(6° • V y ) we have 

div y {(B -Vy^b } = div y {u(a- A)}, div y B = 0. (30) 

Proof, i) For any smooth function u G ker(6°-Vy), taking into account that (& l -V y )j e {o,...,m-i} 
leave invariant ker(6° • V y ), we have 

m— 1 

divj / {(5 • V y u)b } = b° ■ V y (B ■ Vyu) = (V ■ V y u)(b° ■ Vyfa) = div y {u(a - A)}. 

i=0 

Taking now u = ip$ one gets 

b° • V y (B • V^ ) = b° • VyW) = Lb • V y po = L(a - (a )) = (a - A) • V a Vo- 

It remains to establish that LB is divergence free. Notice that for any smooth function 
u G ker(6° • V^), the equation (|29|) becomes 

[b , B] ■ VyU = b • Vy(B • VyU) - B ■ Vy(b • VyU) = (a - A) • VyU. 

In particular taking u G {V'l, V'm-i} we deduce that 

{[B,6 ]+a-A}±span{V^i,...,V^ m _i} = 

implying that 

[B,b°] + a-A = Xb°. (31) 

Taking now u = ipo we obtain 

B • V^ • V v ifo) - b° ■ V y (B ■ Vyi> ) + (a - A) ■ V^o = A/ 

and since we know that 6° • S7 y (B • Vyipo) = (a — A) • VyV'o we deduce that 

B-V y L = XL. (32) 

Taking the divergence of ()3ip implies —6° • Vydiv^i? = 6° • VyA and therefore 

divj / B + A = ^Gker(6°-V y ). (33) 

Combining ([32]), (J33D yields div y (/B) = fJ,L, b° ■ V y n = 0. We are done if we prove that 
divy(LB) has zero average, since, in this case we write 

liL = (/ii) = {divy(LB)) = 0. 

Indeed, by Proposition 13.91 we have 

{m— 1 ^ ( m—1 ^ 

(IB ■ VyiPi) b l \ = div y | < J A> 6*| = 0. 

ii) Obviously we have div^i? = div y B + 6° • VyAo = div y (LB)/ L = and for any smooth 
function u G ker(6° • Vy) we can write 

div y {(i? • V y u)b } = div y {(B ■ V y u)b } = dw y {u{a - A)}. 

□ 
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Remark 3.7 For any smooth function u £ ker(6°-Vj / )U{V , o} w ^ have (B ■ V y u) = (^B ■ V y uJ - 
0. 

Remark 3.8 For any function u € ker(6 • V y ) we have by Proposition \3.1(H 

((a - A) ■ V y u) = (div y (ua)} - (dh y (uA)) = 

and therefore, by Proposition we know that there is a unique zero average function w 
solving b° ■ X7 y w = (a — A) ■ V y u. Proposition \3.11\ says that w = B ■ X7 y u = B ■ V y u. 

Remark 3.9 In our particular case (see Remark \3. 6\) we obtain the coefficients 

qE v v 2 vi qE x v v 

Po = n " Ti2> n = , P2 = — , P3 = n =v ^"n 

mujt: tor uj c uj c muj c \v\ \v\ 



where = ^E/B is the electric cross field drift, and therefore 
B = E = (-— ,»>a) , div x ,«B = - 



V ■ V^c n ( V V " ^x^c _L 

-, B = ,v A k — v 



Q J — I ) "A o 



For any function u(x, v) £ ker T the unique zero average function w(x, v) solving 



is given by 



Tw = (a — A) ■ V x vU = cl ■ V x v u = v ■ V x u H E ■ V v u 

m 



w = B ■ V x v u = B ■ V x v u = j—r ■ V x u + v A ■ V v u. 

L0 C (X) 



The last field we compute will be useful when deriving the second order approximation, in 
Section IP 



Proposition 3.12 Consider a = a(y) a smooth, divergence free vector field in M. m and let B 
be the vector field defined in Proposition ^. Ill ii). Then there is a divergence free vector field 
C = C{y) such that for any smooth function u G ker (6° • Vy) U {ipo} we have 



d\v y {(B ■ Vyu)a}) = div y (uC). 

Proof. For the field B given by (|28p and for any smooth function u G ker(6° • X7 y ) we have 
by Proposition 13.91 

(m l "| 

div y {{B -V y u)a}^ = divjj^ ((B- V y u)(a- V v ^i))6*| 

{m— 1 rra—1 
E E & ■ v v «)e 
i=l j=l 

m— 1 ( m—l ^ 

= E div ^ E u v ■ 
j=i i i=i j 

m—l m—l 

+ ^M^'V^'V.u). (34) 
j=i i=i 
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The last term vanishes, as a contraction between the anti-symmetric matrix ((cti/3j))i t j and 
the symmetric matrix (b l ■ V y (fr J • V y u))ij. Indeed, by Proposition 13.31 we obtain 

= ((on - (ai»&> = ((b° ■ VyPdfr) = - (Pi(b° ■ V y &)> = - (arfi) 

and by Proposition 13.71 one gets 

V • Vyft ■ V y u) -b j ■ V y (6* • V y u) = [6* , V '] ■ V y u = 0, i,j £ {l,...,m- 1}. 

Using again Proposition 13.91 we deduce that for any j S {1, ...,m — 1} 

div y | (oi/Sj) V | = div y | (fta • V^i) &* J = (div v (fta)) . (35) 
Combining ([M]) . ([35]) we obtain 

m— 1 m.— 1 

(divj,{(B • V v «)a}^ = ^ jjV ■ V y u = ^ jjtf ■ V y u 

with the notation 

(div y (/3 + X )Ia)) . n „ ^ w - 

7o = j , 7j = (diVyGSj-o)) , j G {1, ...,m — 1), C=2^7j^- 

i=o 

It is easily seen that the above formula still holds true for u = tpQ. It remains to prove that 
div y C = 0. Since div ?/ (7o6 ) = 0, it is enough to prove that div^j^™^ 1 "fib' 1 } = 0. For any 
smooth function u £ ker(6° • V y ) we have 

— / u div y C dy = C ■ V y u dy = (div y {(.B • V y u)a}} dy = div y {(i? • V y u)a} dy = 
ii m Jm.™ Jr™ ii m 

and therefore (div^C) = 0. But div^C is constant along the flow of 6° 

m— 1 m—1 
div^C = TidiVyfr* + J2b l - V y ji 
i=l i=l 

since 7j,div y 6 J ,6* • are constant along the same flow. Therefore div^C = (div y C) = 0. 
□ 

Remark 3.10 Actually we will need a divergence free vector field C satisfying the statement 
in Proposition \3.1S\ only for functions u 6 ker(6° • V y ). The coordinates along (6 J ) ig { li ... im _n 
of such vector fields are uniquely determined by 7, = (divy(^a)), i £ {l,...,m — 1} and 
the coordinate along b° could be any smooth function in ker(6° • Vy), which guarantees the 
constraint div y C = 0. 

Remark 3.11 In the particular case of the vector fields b° = (0, 0,uj c (x) v), 6 1 (1, 0,0,0), 
b 2 = (0,1,0,0), b 3 = (0,0,v/\v\), a = (v,q/mE) we obtain the coefficients 

f x \v\va V X B 
(71,72) = -v A - v GD , 73 = - — 

where Uaj^GD are the (rescaled) electric cross field drift and magnetic gradient drift 

± E \v\ 2 L V X B 



B ' ^ u 2w c (x) B 
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and the vector field C becomes 

C = (-v A - vgd, - VA 2^ v + ^o^cix) ± v) (36) 

for some 70 6 kerT. It is easily seen that div x u A = — (u A • V X B)/B, div x fGD = and thus 
div XjV C = 0. In the sequel we will consider 70 = 0. 

Proposition 3.13 Consider a = a(y) a smooth, divergence free vector field in W m . Then for 
any smooth function u we have 

div y {(B ■ V y u)A} = div v {u[A, B]} + div y {(A ■ V y u)B} 

where the vector fields A and B are given by ([27)) and Proposition VJ.11\ ii), respectively. 

Proof. The field [A, B] is divergence free 

div y [A, B] = A ■ Vydiv y B - B • V y div y A = 

and therefore 

div y {(B ■ V y u)A} = A-V y {B -V y u) 

= [A, B] ■ V y u + B ■ V y (A ■ V y u) 

= div y {u[A,B]} + div y {(A-Vyu)B}. 

□ 

Remark 3.12 Let u G ker(6° • V^) and z satisfying 

div y {u(a - A)} + dw y {zb } = 0, (z) = 

which means that z = —B ■ VyU, cf. Proposition \3.11l Then combining Propositions \3.12\ 
and \3.l3\ we obtain 



{divy(za)) - d'wy(zA) = - (^divy{(B ■ V y u)a}J + div y {(B ■ V y u)A} 

= -div y (uC) + div y {n[^, B}} + div y {(A ■ V y u)B} 
= dw y {u([A,B]-C) + (A-Vyu)B}. 

3.3 Average and the collision operator 

In this section we will show that the average operator associated to T = oo c (x) ± v ■ V„ and 
the collision operator are commuting. This means that at the lowest order the collision mech- 
anism is not affected by strong magnetic fields. Corrections appear only at the next orders, 
as commutators between transport and convolution through the scattering cross section, see 
Proposition 13.151 We first recall the classical properties of the linear Boltzmann operator Q. 
The notation L 2 M stands for L 2 . 

Proposition 3.14 Assume that the scattering cross section satisfies ©. Then the gain/loss 
collision operators Q± 

Q+(f)(x,v) = - f s(v,v')M(v)f(x,v')dv', Q.(f)(x,v) = - [ s(v,v')M(v')f(x,v) dv' 



are linear bounded operators from I? M to I? M and ||Q±|| < Sq/t. Moreover, if the scattering 
cross section also satisfies (jU), then Q± are symmetric with respect to the scalar product of 



T 2 
M' 
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Notice that M depends only on \v\ and therefore M G kerT, T = oj c {x) v ■ V v . Observe 
that the average operator is also bounded from L 2 M to l? M . Indeed, for any / G I? M we have 
f/y/M G L 2 (dvdx) and thus 



|2 



f 



< 



L 2 



f 



L 2 



2 



We assume that there is a function a : R + — > [so, So] such that 

s(y,v') = a(\v - v'\), v,v' G R. 



(37) 



Lemma 3.1 Assume that the scattering cross section satisfies ()37p . Then the gain/loss 
collision operators are commuting with the characteristic flow ofT. In particular the gain/loss 
collision operators leave invariant L M n kerT and L 2 M n ker (•). 

Proof. Consider / G l? M and for any a G R let us introduce f a (x,v) = f(x,R a v), where 
R(a) stands for the rotation of angle a (see (fI7))). Obviously / Q G and 

(Q+ZaJC^u) = -/ a(\v-v'\)M(v)f a (x,v')dv' 



= (Q+f)(x,R a v). 

Similarly one gets 

(Q-f a )(x,v) = -f a(\v-v'\)M(v')f a (x,v)dv' 



-/ a(\v - v'\)M{v)f(x,R a v') dv' 

<r(\v - R~ a w'\)M(v)f{x,w') dw' 
a(\R a v -w'\)M(R a v)f(x,w') dw' 



-I a{\v -v'\)M(v')f(x,R a v) dv' 

a(\R a v - R a v'\)M(R a v')f(x, R a v) dv' 



a(\R a v - w'\)M(w')f(x, R a v) dw' 



= (Q-f)(x,R a v). 

In particular, if / G L 2 M n ker T then 

(Q±f)(x,R a v) = (Q±f a )(x,v) = (Q±f)(x,v) 

saying that Q±f G kerT. Consider now h G Linker (•) and ip G kerT Then Q±(h) f\~M G 
L 2 (dfdx), y/~M(p G L 2 M and we can write 



Q±(h) 



<pdvdx = (Q±(h),VM<p) L 2 u = (h,Q±(VM<p)) L 3 



h Q±(VMp) 



dvdx. 



L 2 n ker T Since h/\/M G L fl ker (•) we deduce that 



By the previous computations we know that Q±{\fMip) G Linker T and thus Q±(y/~M<p) / \/~M G 



dudx = 



□ 



Q±(h) 



h Q±(VM<p) 



03 dvdx - 
'M Jr 2 Jr 2 VM x/M 

saying that ( Q±(h)/VM) = and thus Q±(h) G L 2 M n ker (• 
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Corollary 3.3 The average and the gain/loss collision operators are commuting on L 2 M . 

Proof. Consider / G L? M . Using the decomposition / = (/) + (/ — (/)) one gets by the 
previous lemma that TQ± (f) = 0, (Q±(f — (/))) = and therefore 

<Q±(/)) = (Q± (/)) + (Q±(f - (/))) = Q± (/) • 

□ 

We end this section with the following commutation result, which will be used when investi- 
gating the second order approximation. 



Proposition 3.15 Assume that the scattering cross section satisfies (|37p and it is smooth. 
Then for any smooth function f = f(x,v) £ kerT" we have 



Q\f) := Q(B-V XjV f)-B-V X)V Q(f) 

= ~ [ (si- tjs 2 ) M{v)—^ ■ V x f(x, v') dv' 

m , s f , . h / , , m , . f It/ 1 



t0 J R 2 t6 J r2 \v\ 



where (see Remark \37l 

~ ( ± v v-V x oo c \ . , ... . / . v ■ v' \ 

B=[ ,v A = — v), si(v,v) = (s(v,v)), s 2 {v,v ) = { s{v,v ) - ) ■ 

\ UJ C lj£ J \ \v\ \v'\ I 

Proof. Observe that V^s + V„/s = and therefore, performing integration by parts with 
respect to v' yields 



tB ■ V x , v Q(f) = / s(v, v'){M(v)B(x, v) ■ V x yf(x, v') - M(v')B(x, v) ■ V x , v f(x, v)} dv' 
Jr 2 

+ B v (x,v)- [ {V v M(v)f(x,v')-V v/ M(v')f(x,v)}s(v,v')dv'. 

JR 2 

Combining with the equality 

tQ(B ■ V x ,«/) = f s(v, v'){M(v)B(x, v') ■ V xy f(x, v') - M(v')B(x, v) ■ V x>v f(x, v)} dv' 
Jr. 2 

we deduce that 

rB-V XtV Q(f) - rQ(B-V x>v f)= [ sM(v)(B(x, v) - B(x, v')) ■ V x yf(x, v') dv' 

Jr. 2 

+ B v (x,v)- I {V v M(v)f(x,v') -V v ,M(v')f(x,v)}s dv'. (38) 
Jr 2 

We transform the above terms appealing to the symmetry of /. Indeed, we have 



s(v, v')M(v) ^ . • V x f{x,v') dv' 



s(v,v')M(v)(B x (x,v) - B x {x,v')) ■ V x f(x,v') dv' 
{v'-t 

U> c (x) 

^4 ( s ( v > v ') V ~ v )) ■ dv'. (39) 

:U C {X) \ I 
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We introduce the following averaged (with respect to v') scattering cross sections 

si(v,v') = (s(v,v')) , s 2 (v,v') = (s(v,v') j^j~~j^7j^ • 
Writing that v' = r'cosa v/\v\ — r'sina _L t'/|f|,r / = \v'\,r = \v\ one gets 



v — v = (r cos a — r) 



V , V 

— — — r sm a -r—r 

\v\ \v\ 



1 



2~ 



si(v,v') = — / c(v (r 1 ) 2 + r 2 — 2rr' cos a) da 
2vr Jo 

S2(v,v') = — / c(v (r 1 ) 2 + r 2 — 2rr' cos a) cos a da 
27r Jo 

saying that s±,S2 depend only on r = \v\,r' = \v'\ and that they are symmetric with respect 
to r,r'. Notice also that 



I r 2n 

(s(v,v')v') = — / o"(v ( r ') 2 + r2 ~~ 2rr'cosa) 



r cos a-r-r — V sm a-p-r 



da = S2\v'\ 



(s(v, v')(v - v')) = ( si(v, v') - — s 2 (V) v') j v. 



and 

Therefore, the integral on the last line of (|39p writes 

Ti = - f (si(v, v') - ^-s 2 (v, v')) M(v)— ^ ■ V x f(x, v') dv'. 
By taking into account that v' ■ V v if(x, v') = 0, we observe that 

T 2 := [ s(v,v')M(v)(B v (x,v)-B v (x,v'))-V v ,f(x,v')dv' 
s(v,v')M(v) 



(40) 



s(v, v')M{v) V ' V f c d\v v ,{ L {v' - v)f(x, v')} dv' 



-' ' 1 " ' ■•' 1 ' " ■ dv' = 0. 



\v — V i 

For the last integral in (|38p we notice that 



A\ v -^\)jl r ± (v'-v)f(x,v')M(v)- 



(41) 



T 3 := B v (x,v)- I V v M(v)f(x,v')s(v,v') dv' = ~(v-v A ) I Sl (v,v')M(v)f(x,v') dv' (42) 



and 



-B v (x,v)- [ V v/ M(v')f(x,v)s(v,v' ) dv' 
Bv(x,v)- / -(s(v,v')v')M(v')f(x,v)dv' 



m (vA-v) I S2(v,v') l -^-M(v')f(x,v)dv'. 



(43) 
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Finally combining (J38|), (00]), (JHD, g2]), (g3]) we obtain 

OH/) = -/ (si(v,v>) - l -fls 2 (v,v')) M(v)v ■ \^v A f(x,v') - ±V *fM ) dv' 
t Jr2 V \v\ ) [9 ui c {x) J 

I f I V I Tfi 

+ - s 2 (v,v'y—l-(vv A ){M(v)f(x,v , )-M(v')f(x,v)}dv' 
T JR2 \v\ 9 

= If ( Si ^, v ' ) -\1L1 S2{VjV ' ) )m(v)F(x,v')^-V x (§^) dv' 
T J R 2\ \v\ J uj c {x) \F(x,v')J 

t Jrz \v\ 9 { 



F(x,v') F(x,v) 



dv' 



where F(x,v) = M(v)ex.p(—q(j)(x)/9) is the equilibrium of ([7]), E(x) = —V x (j). \ — | 
Notice that the operator Q 1 has the following properties 

1. the equilibrium F belongs to the kernel of Q , since B ■ V XV F = and Q(F) = 

QHF) = Q(B • V XjV F) - B ■ V X , V Q(F) = 0. 

2. Q 1 maps constant functions along the flow of b° to zero average functions since, for 
any smooth function / G kerT, B • V x>v f £ ker (•), Q{f) G kerT cf. Lemma I3TT1 and 
therefore 

Q(B ■ V x ,vf) e ker (•) , B ■ V x , v Q(f) G ker (•) . 

3. In particular Q 1 satisfies the local mass conservation i.e., for any smooth function 
/ G kerT we have Q l {f) G ker (•) and thus 



Hf)dv= [ (Q\f)) dv = 0. 



4. Q 1 is symmetric with respect to the scalar product of L F = L 2 y p^ X v ) J • Indeed, for 
any smooth f,g£ ker T we have 

(QHf),9)Ll = I [ Q(B-V x , v f)^dvdx- [ I B-V X , V QUA dvdx 

= (Q(B ■ V x ,vf),g)L% + (Q(f),B ■ V x , v g) L% 
= (B ■ V x , v f, Q(g)) L% + (Q(f),B ■ V x , v g) L% 
= {Q\g),f) L 2 F . 

5. For any smooth / G ker T we have (Q (/), /)l 2 = since Q 1 (f) G ker (•), f/F G ker T. 

Remark 3.13 When the scattering cross section is constant, we have si = s,S2 = and the 
expression of Q 1 becomes 

Q\f) = S -M(v) PF {x)^- ■ V x (El) 

T UJ C (X) \PF J 

where p f (x) = J R2 f(x,v') dv' , p F (x) = j R2 F(x,v') dv' = exp(-#(x)/0). 
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4 Asymptotic models for Boltzmann equation 

This section is devoted to the derivation of the limit model for d?D when e becomes 
very small. At the formal level, we perform our analysis starting from the expansion f e = 
f + ef 1 + e 2 f 2 + .... We are interested in first order but also second order asymptotic models. 
We present weak and strong convergence results. 



4.1 First order model 

We investigate the asymptotic behaviour of the linear Boltzmann problem ([7J, ([2]) when 
e \ 0. At least formally the leading order term / satisfies ([9]) and the constraint Tf(t) = 
0, t E R+. By Proposition 13.21 we know that Range T = ker (•) and therefore © is equivalent 
to 

(d t f + v-V x f + ^-E- X7 v f) = (Q(f)) . 
\ ml 

It is easily seen that {dtf) = dt (/) = dtf and by Corollary 13.31 we know that (Q(f)} = 

Q((f)) = Q(f)- Moreover, there is a function g such that f(t,x,v) = g(t,x,r = \v\) and 

therefore 

(v ■ V x f) = (v) ■ V x g = 0, (E- V„/> = (e ■ d r g\ = yj ■ E d r g = 

\ \ v \ I \ v \ 

since {v ) = 0. We obtain the limit model 

dtf = QUI Tf = 

which can be justified rigorously at least under the confinement potential hypothesis 
E(x) = -V x (j){x), [ exp f-q^l] dx < +oo. 



Indeed, in this case F(x,v) = M(v) exp (— qcp(x)/9) belongs to 1? F and solves ([7]) for any 
e > 0. Notice that Q± are also bounded from l? F to Lp. It is well know that for any convex 
function H, the relative entropy 

H E (t)= [ [ F(x,v)H ( f y/ X,V A dvdx 



F(x, v) 
satisfies 

-jL+V s (t) = 0, t G M + (44) 
where the dissipative term T> 6 is given by 

V e {t) = -( f [ s(v,v')M(v)F(x,v')D ( f y/ XlV x \ f y/ X, f^ ) dv'dvdx (45) 



t Jr 2 Jr 2 Jr 2 ' ' V F(x,v) ' F(x, v') 

where 

D(u, v!) = H{u') - H{u) - (u' - u)H'{u) > 0, u, u' G M. 
Assuming that the initial condition satisfies 

r>o,re4 (46) 

one gets by (|4"4"]) 



F(x,v) ' ' Jr2J R 2 F(x,v 
which allows us to pass to the limit, at least weakly, in ([7|). 
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Proof, (of Theorem I2.ip We have the uniform bound 

SUP 1 1 f |Loc (M L 2) < ||/ in || L 2 

£>0 V F ' F 

and we can assume (after extraction of a sequence \ 0) 

lim f £k = f, weakly * in L°° (R+; L 2 F ). 

Using the weak formulation (J7J) it is easily seen that the weak limit / satisfies the constraint 
Tf = 0. Taking now test functions of the form T](t)<p(x,v), tj G C£(R+), ip G C£(R 2 x R 2 ) n 
ker7~ (the constraint cp G ker T is imposed in order to remove the singular term \Tf e in (JJJ)) 
we obtain 

- / / / f k (v'<P + rjv ■ V x <p + r]—E ■ V v ip) dvdxdt - /N(0)v dvda; ( 47 ) 

Jr+Jr 2 Jr 2 v m -' Jr 2 Jr 2 

Q(/ £fe )w dvdxdt. 

2 

By the symmetry of the collision operator one gets, noticing that rjQ(Fip) G L X (R + ; L 2 ^) 
lim / / / Q(f*) W dvdxdt = lim / (/ e *(t), V (t)Q(F<p)U . di 

(f(t),r,(t)Q(F<P))i%M 
Q(f)rj<f dvdxdt. 

2 

Passing to the limit with respect to — > +oo in (|47p leads to the weak formulation 

- / / / / (v'f + Vv ■ V x tp + V— E ■ V v <p) dvdxdt - f ' m v(Q)<P d vdx (48) 

Jr+JbizJe. 2 v m J Jr2Jr2 



Q(f)V<P dvdxdt 

2 

for q G C£(K+)> 92 G C^R 2 x M 2 ) n kerT. Writing tp(x,v) = ip(x,\v\) we deduce that 
(v ■ S7 x <p) = 0, (E ■ V v ip) = and therefore, since Tf = 0, one gets 



f [rjv ■ V x tp + rj — E ■ V v <p dvdxdt = 0. 
2 V m 

Finally the formulation (|48|) reduces to 

fdtirup) dvdxdt - f [ </ in > 77(0)^ dvdx = III Q{f)w dvdxdt 



for test functions satisfying the constraint T(j]ip) = at any time t G R+. Since /, (/ m ) 
and Q(/) also satisfy this constraint it is easily seen that / solves weakly (j 1 3 H . By the 
uniqueness of the solution for (fT5j) we deduce that all the family (f £ ) £ >o converges weakly * 
in L°°(M_|_; L F ) towards /. It remains to justify that the constraint Tf(t) = is automatically 
satisfied at any time t > 0. Indeed, for any a > we have, cf. Lemma 13.11 

dtf a = (d t f)(x,R a v) = (Qf)(x,R a v) = (Qf a )(x,v). 

Since f a (0,x,v) = (f m ) (x,R a v) = (/ m ) (x,v) = /(0, x,v) we deduce that both / and f a 
solve (flUj) and, by the uniqueness of the solution, one gets f a (t) = f(t) for any a G R, saying 
that Tf(t) = at any time t > 0. q 
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Under appropriate smoothness hypotheses on the initial condition f m and the electro- 
magnetic field, it is possible to obtain strong convergence results. Actually we can show that 
f £ = f + 0(e) strongly, which justifies the first term in the expansion f £ = f + ef 1 + ... . We 
only sketch the main lines of the arguments and skip the standard details. The main idea is 
to introduce a first order correction term. Assume for the moment that the dominant term 
/ is smooth enough (with respect to all variables). Since 

(d t f + v-V x f + ^-E- V v f - Q(f)\ = 
\ m I 

we deduce by Proposition 13.21 that there is a unique function f 1 such that at any time t > 

d t f + v-V x f + ±E-V v f + Tf 1 = Q(J), </ 1 > = 0. (49) 
m 

Combining with ([7]) we obtain 

( d t + v ■ V x + ±E ■ V v + -) (f-f- ef 1 ) = Q(f - f - ef 1 ) + eKtf 1 ) (50) 
\ m e J 

where 

W 1 ) = -(dt + vV x + ±E- V v ) f 1 + Qif 1 ). (51) 
V m / 

Notice that lZ(f l ) is of order of 1, provided that Z 1 is smooth enough (we will see later on 
how we can get regularity for f 1 from the regularity of /). Assume for the moment that f 1 
is smooth enough. In this case e1Z(f 1 ) will be of order of e. Using again the relative entropy 
dissipation (similar to (|44|) . since now we have an extra term elZ(f 1 )) or multiplying (|5U|) by 
(f e — f — ef 1 )/F, where F is the equilibrium of (JJJ), one gets 

< e\\KU\t))\\ L *J(r - f - ef l )(t)\\ Ll 



implying that 



Finally one gets 



^ii(r-/-e/ 1 )(*)ii^<«/ i (t))n^ 



ik/ £ - / - efmhi < nr - </ in ) - efrnki +ef t \\nf\s))\\ Ll d S 

and therefore, for any T > 0, there is a constant Ct depending on ||/ 1 ||Loo( K+;i |) and 

IIW 1 )ll£i(0,T;X|,) such that 

\\(f £ -f)(t)\\ L2 <C T e, te[0,T], e>0 

provided that the initial condition is well prepared i.e., Tf m = 0. It remains to justify the 
smoothness of / and / . 



Proposition 4.1 Assume that the initial condition satisfies f m , \v \f m G L F , Tf m = 0. Then 

-5 Lp). 



the solution of $T3D satisfies f G L°° (R+; L 2 F ) , vf G L°°(lR + ;Lf 
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Proof. We know that ||/(i)||x,2 < ||/ m ||L|,- Multiplying the equation in (fl"3|) by v one gets 

dt(vf) + vf(t,x,v)- [ s(v,v')M(v') dv' = vM(v)- [ s(v,v')f(t,x,v')dv'. 

Solving the above differential equation with respect to vf and estimating the right hand side 
by using 

/ s(v,v')f(t,x,v')dv' <S \\r\\ L 2 

lead to the bound 

\\vf(t)\\ L2 <\\vn\ L2 I \vfM(v')dv') 1/2 \\n\ L2 *>o. 

F F SQ \J R 2 J F 

□ 

Proposition 4.2 Assume that the function a, introduced in (|37p . is Lipschitz continuous 
and that the initial condition satisfies 

rr = 0, (i + \v\){\r\ + \v x r\ + \v v r\) G l\. 

Then the solution of (TT3"j) satisfies 

(1 + |«|)(|/| + \d t f\ + \V x f\ + \V v f\) G L°°(1R + ; L 2 F ). 

Moreover, if a G W 2 >°° and (1 + M 2 )V^/ in e L| ifcen (1 + \v\ 2 )V 2 t ^ v f G L°°(R+ ;!$.)■ 

Proof. By Proposition 14.11 we know that (1 + \v\)f G L°°(M.+ ;Lp) and since the operators 
Q± remain bounded in L 2 F (see Proposition 13. 14j) 

9 Q 

||Q(/)IL| < — 

F f F 



one gets 11^/11^(^.^2) < ^H/" 1 ^ . Observe that 

s(v,v')f(t,x,v')dv'< ( [ s(v,v') 2 M(v')dv') 1/2 \\f(t,x,-)\\ Ll <S \\r(x 



Therefore we obtain 

1 /2 

\\vd t f\\ L 2 p = \\vQ{f)\\ L 2 F < ^ \^Jjv'\ 2 M{v') dA \\r\\ LF + 

saying that vdtf E L°°(IR + ; L 2 ^). Since Q is commuting with the space derivatives one gets 

dt(V x f) = Q(V x f), (t, 1 ,«)ei;xl 2 xl 2 , V^OHV^GkerT 

and therefore, by Proposition 14.11 we obtain (1 + H^V^/I G L°° (R+; L 2 F ) . It remains to 
estimate the velocity derivatives. As we already know that • V„/ = let us analyze the 
derivative along v. We have 

dt (^t • V v f) +-(^r- V v f) [ s(v, v')M(v') dv' = -g (52) 

\\V\ J T \\v\ ) J R2 T 
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where 



,; / ■ V v s{M(v)f(t, x, v') - M(v')f(t, x, v)} dv' + ^j- V v M(v) \ s(v, v')f{t, x, v') dv' 
/r 2 M 



Observe that 



g\\i><2\\o'\\ L ~\\r\\iz+s \\r\\iA 



I 2 

M 



and finally, solving ([52]) with respect to A • V v f yields 



V„/(i) 



< 



; 2 

F 



H [ 2||ct / ||loo + 5o 

Ll S \ 



V 




—r ■ V V M 


J 


\v\ 


M J 



wrhi 



Therefore V v f G L°° (R + ; L 2 F ) . Starting from the equality d t (v ■ V v f) = v ■ V v Q(f) we obtain 
the differential equation 

dt(v ■ V v f) + -{v ■ V v f) [ s(v, v')M{v') dv' = -h 
t Jr2 r 

where 

h= f v V v s{M(v)f(t,x,v') - M(v')f(t,x,v)} dv' + v- V V M [ s(v,v')f(t,x,v')dv'. 



Observing that 

Mt)\\ Ljp < llo^lU.-llvArili^ll^llx,- -i- II^IU~ll«.fC*)llx,» H-^sbll^ - v^Aril^^n^n^a 

we deduce that \\v-V v f(t)\\ L 2 F < \\v-V v f m \\ L 2 F +1/so||^||i,«'(r + ;L 2 ) which says that \v\ \V v f \ = 
\ v " V v /| G L°°(lR + ;Lp). The last statement follows similarly. |— I 

Proposition 4.3 Assume that the hypotheses of Propositions EOl and \4~J^ hold true and that 
the electro-magnetic field (E(x),B(x)) belongs to W 1,co (M. 2 ). Then we have 

f l eL°°(R + ;L 2 F ), 1Z(f l ) G L°°(R + ;L F ) 

where f 1 is defined by (gSJ and TZif 1 ) by ([51]) . 

Proof. Since dtf = Q(f), notice that (09]) can be written 

v ■ V x f + ^-E(x) ■ V v f + r/ 1 = o, </ 1 > = o 
m 

and therefore we have, cf. Remark 13.91 
f 1 = -B ■ V* ,„/ 
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w c (a;) raw c (i) 
Our conclusions follow easily by the smoothness of /. 

Combining the previous arguments we have proved the strong convergence result stated in 
Theorem 
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4.2 Second order model 



As we have seen in the previous section, the leading order term in the expansion f £ = 
f + ef 1 + ... satisfies the space homogeneous Boltzmann equation. The collision operator 
appearing in the first order limit model is still the linear Boltzmann operator. In other words 
the magnetic field does not affect the collision mechanism (the collision operator and the 
advection along the Laplace force oj^tv ■ V„ are commuting). Our goal in this section is to 
take into account the first order correction. We are looking for a closure involving / + ef 1 
accounting for the small departures from the equilibrium / when the magnetic field grows 
like 0(l/e). We expect that the magnetic field will affect the collision mechanism. We will 
see that the new collision operator entering the second order model is no more local in space. 
Naturally, the new collision operator should verify the usual physical properties like mass 
conservation, entropy dissipation. The derivation is completely formal and it is based on the 
average technique introduced before. The limit model being quite complicated we prefer to 
skip the technical details leading eventually to rigorous results. Nevertheless we expect that 
the arguments employed when studying the first order model will adapt to the second order 
one as well. 

The linear Boltzmann equation (J7|) reads 

dtf + a ■ V X)V f + -Tf = Q(f ) (53) 

where a = (v,q/m E(x)),T = uj c (x) ± v • V«, w c (x) = qB{x)/m. Plugging the expansion 
f £ = f + ef 1 + ... in (f53j) yields the constraint Tf = and the evolution equations 

dtf + a- V XiV f + Tf 1 = Q(f) (54) 

d t f l +a-V x , v f l + Tf = Q(f l ) (55) 



The evolution equation for the dominant term / has been obtained by averaging (|54p along 
the flow of T ■ Indeed, by Proposition 13.101 we know that there is a vector field A = A{x, v) 
such that 

(a • V x ,„/) = A • V^/, / GkerT. 

Actually the average field associated to a = (v,q/m E(x)) is A = cf. Remark 13.61 but 
we prefer to keep track of the vector field A, which will clearly show that our technique will 
provide a second order closure independently the field A vanishes or not. Nevertheless, in 
the final model the field A will be replaced by its value ^4 = 0. For any / S ker T we have 
(dtf) = dt (/) = dtf and by Corollary 13.31 we can write 

<Q(/)> = Q«/» = Q(f). 

Therefore averaging (|54p leads to the homogeneous Boltzmann equation 

dtf + A-V XtV f = Q{f). (56) 

Clearly more information is available from (|54f) . In order to exploit that we use the decom- 
position 

f^g^h 1 , g 1 =(f) GkerT, h 1 = f 1 - (f 1 ) G ker (•) 
and therefore (|54p becomes 

dtf + a-V^f + Th 1 =Q(f). (57) 
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The equations ([56]) . (f57|) yields 

(a- A)- V x , v f + Th l =0 (58) 
and by Remark 13.91 we deduce that h 1 = —B ■ V x ,vf where the field B writes (see Remark 

B = — ^- -v, — L E v 



uj c {x) \ ' m w c (x) 
We average now (|55|) . Notice that we have 

(d t f 1 ) = d t (f 1 ) = d t g\ (Tf 2 ) = 0, {Q{f 1 )) = Q{{f 1 )) = Q(g 1 ) 

(a • V^/ 1 ) = (a • V^g 1 ) + (a • V^/i 1 } = .4 • V^,^ 1 - (a • V^(B • V^/)j 
and therefore we obtain the equation 

V + .4 • V^g 1 - (a • V X>V (B ■ V,,J)) = Q^ 1 ). (59) 

Naturally, the closure for / + ef 1 = f + eg 1 + eh 1 comes by combining (f56]h ([59]) 

+ 4 • V x>v )(f + e/ 1 ) + s((o • V I)B ) - 4 • V X:V + Q- d t )h l = Q(f + ef 1 ). (60) 

Making use of Remark 13.121 and Proposition 13.151 we obtain 

((a • V*,„) - 4 • V x , v + Q~ d t )h l = -div x , v (fC) - Q(B ■ V x>v f) + d t (B • V x , v f) 
= -div x>t ,(/C) + B • V,,«Q(/) - Q(B • V x ,«/) 
= -C-V x , v f-Q\f). 



Finally (|60|) becomes 

+ (4 - eC) • V,,„)(/ + e/ 1 ) = Q(f + e/ 1 ) + eQ\f) - e 2 C ■ V^/ 1 (61) 

and since we are looking for a closure with respect to / + sf , or something equivalent 
to / + e/ 1 up to a second order term 0(e 2 ), we need to define Q 1 over the space of all 
functions (not only on the subspace of constant functions along the flow of 6°, as it was done 
in Proposition 13.15]) . We denote by Q 1 this extension and assume for the moment that Q 1 
is known (the exact expression will be determined in Proposition 14.41 and Corollary 14. 1|) . By 
doing that we observe that the right hand side of (|6ip writes 

(Q + eQ l )(f + ef 1 ) - e 2 (C ■ V^f 1 + Q 1 (f 1 )). (62) 

Getting rid of the second order term and taking into account that A = 0, we may expect 
that solving 

d t f - eC ■ V x , v f = (Q + eQ^Cf) (63) 
provides a second order approximation of f £ , since at least formally 

f ~f e = (f - (/ + ef 1 )) - (e 2 f 2 + -..) = C(e 2 ). 

A delicate point is how to pick the extension of Q . Naturally, the global collision operator 
Q + eQ 1 , and also the first order correction Q 1 , should dissipate the entropy 

(QHf)J)Li<o, v/. 
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In particular, for any smooth functions g G ker T, h G ker (•) and any 5 > we have 

(Q\ g + 5h),g + 5h) L , F <0. (64) 

Recall that one of the properties of Q 1 was that (Q (#), 5) r2 = and therefore becomes 

(Q 1 ^),^! + (g,Q\h)) L% +8(Q 1 (h),h) L 2 F < 

which is equivalent to 

{Q l (g),h) L 2 F +(g,Q 1 (h)) L 2 F = 0, (Q 1 ^), %2 < 0, g G ker T, /» G ker (•) . (65) 

Recall also that Q 1 maps constant functions along the flow of b° to zero average functions. 
We will construct the extension of Q 1 such that Q 1 maps zero average functions to constant 
functions along the flow of 6°. Notice that in this case the inequality (Q l {h), h) L 2 < holds 
true for any h G ker (•) (actually it becomes equality, by the orthogonality between ker T and 
ker (•)) and that the first condition in (|65p uniquely defines Q 1 (/i) G ker T for any h G ker (•). 
Moreover, the above construction will ensure the global mass conservation since Q 1 satisfies 
the local mass conservation and for any smooth zero average function h we can write, recalling 
that Q l (F) = 

/ / Q 1 (h)(x,v) dvdx = {Q\h),F) L 2 = -(h,Q\F)) L 2 = 0. 



Proposition 4.4 For any smooth functions g G ker T, h G ker (•) we have 

1 



(QHg),h) 



i 2 

F 



i 2 

F 



where 



K{h) 



Jr*\ \v'\ J 



h(x, v') 
u c (x) 



dv' 



m 
s 2^r 



\v\M(v)^—^ h(x,v') - \ v '\M(v')^-^ h(x,v) 



dv'. 



In particular the unique function Q l {h) G ker T satisfying 

(Q 1 (g),h) Ll + (g,Q 1 (h)) L 2=0, V<?GkerT 



is given by 
Q\h) 



si 



jrs 2 I M(v)d\\ x 



~v'h) (x, v') 

LO c (x) 



dv' 



m 



S2- 



. (v'h)(x,v') . „,,. (vh)(x,v) 
\v\M(v) v A ■ - , ' ; - \v'\M(v') v A ■ * — 1 \ - 1 



(66) 



dv'. 



Proof. By Proposition 13,151 we know that 

Q\g) = - 



W\ \ , . ,x v „ / q(x : v') 



+ -/ rjS 2 -(vv A )M(v)F(x,v') 

T Jr 2 \V\ 



g{x,v') _ g(x,v) 
F(x,v') F(x,v) 



dv' 
dv'. 
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Integrating by parts with respect to x leads to 



r(Q\g),h) Li 



v 



h(x - V ^ Six -" r > dv'dvdx 



It? I Tfl hi x 

Ir?Jr? \v\ F(x,v) 



-(g,K(h)) L 2 F . 



Since 



r(g,Q\h))v F = -T(Q\g),h) L 2 F = (g,K(h)) Lp = (g,(K(h))) Lp 

and Q l (h) S ker7~, we deduce that Q l {h) = (K(h)) jr. For further transformations of the 
terms in K(h) we need the easy lemma 

Lemma 4.1 Let k = k(v,v') be a bounded function depending only on r = \v\ and r' = 
\v'\. Then for any function f = f(x,v) S l? F we have J R2 k(v,v')M(v)f(x,v') dv' G L 2 F , 
J R2 k(v , v')M (v') f (x, v) dv' £ l? F and 

k(v,v')M(v)f(x,v') dv') = [ k(v,v')M(v)(f)(x,v')dv' 



k(v,v')M(v')f(x,v) dv' 



k(v,v')M(v') (f) (x,v) dv'. 



Proof. Observe that J R2 k(v,v')M(v) (f) (x,v') dv', J R2 k(v,v')M(v') (f) (x,v) dv' depend 
only on x and \v\ and therefore they belong to kerT. Pick now (p = (p(x,v) £ L F a constant 
function along the flow of 6°. We can write, observing that (x, v') — > J R2 k(v, v')M(v)(p(x, v) dv S 
ker T 



[ <p [ k(v,v')M(v)f(x,v') dv' 
Jr 4 Jr 2 



dv dx 




f(x,v') / k(v, v')M(v)ip(x, v) dv dv'dx 
2 Jr 2 

(/) (x, v') I k(v, v')M(v)tp(x, v) dv dv'dx 
2 Jr 2 

(p(x,v) / k(v,v')M(v) (f) (x,v) dv' dvdx 



saying that (f R2 k(v,v')M(v)f(x,v') dv') = f R2 k(v,v')M(v) (f) (x,v') dv'. The second state- 
ment follows immediately since f R2 k(v,v')M(v') dv' 6 kerT and thus 

/ k(v,v')M(v')f(x,v) dv') = [ k(v,v')M(v')dv'(f)(x,v). 
Jr 2 I Jr 2 



□ 



Appealing to Lemma 14.11 and Remark 13.41 we obtain 



and 



m 



m 



\v\M[v) (-^ • v A J h(x, v') - \v'\M(v') (-^ ■ v A J h(x, v) 
\v\M(v)v A ■ ^Oj^/ i x ' v ') ~ \v'\M(v')v A ■ (rih) ( x i v ) 



dv' 
dv' 
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implying that 



Q\h) = - [ ( Sl - ^rs 2 ) M(v)div 



V/i)(x,u')\ , 

X I — I CLV 



u c (x) 



III 



•52" 



, _,. s (v'h)(x,v') . /.,,- , N (vh) (x,v) 
\v\M(v) v A ■ - , ' ; - \v'\M(v') v A ■ - — , ' 1 



dv'. 



□ 



Corollary 4.1 The extension Q 1 of the operator Q 1 satisfying (J66J) is given by 

Q (/ = - / si - yjS2 M(v) v ■ —— {f) (x,v) n dv 

(v ■ v A ){M(v) (/) (x,v r ) - M{v') (/) (x,v)} dv' 



If \v I m 
+ - / S 2 y^ 



1 

+ - 
T 

1 



\v'\ J \ U) C (X) I 



m 

S2-z v A ■ 



\v\M(v) {V ' f ]^ V,) - \v'\M(v>)^^ 



dv'. 



Moreover, when the scattering cross section is constant we have 



Q\f) = -M(v) 



div-r 



uj c (x) ) lo c (x) x \pf{x) 



+ 



(67) 



(68) 



where p f = J R2 / dv', p F = J R2 F dv' , j f = J R2 v' f dv' . 

Proof. Using the decomposition / = (/) + (/ — (/)), with (/) 6 ker T and / — (/) G ker (•) 
and taking into account that (v(f — (/))) = (vf) we can write 

Q\f) = Q 1 ((f)) + Q\f-(f)) 

leading to the formula (|67p . In particular, when the scattering cross section is constant, we 
have s\ = s, S2 = and ([67]) reduces to (j68|) , | — | 

Remark 4.1 Replacing the vector field C by (|36p . the collision operator Q by (J3l) and the 
extension operator Q 1 by (j68j) in the equation (|63p leads to the second order model stated in 
Theorem \2.3\ (case of constant scattering cross section). 

It is easily seen, at least for constant scattering cross section, that the kernel of Q + eQ 1 
reduces to MF = {XF : A G M}. Let us establish this property in the general case. 

Proposition 4.5 For any e > the only functions in the kernel of Q + eQ 1 are those 
proportional to F. 

Proof. Let / be such that (Q + eQ 1 )(f) = 0. Based on the properties of Q and Q 1 we have 

= (Q(f) + eQ\f)J) L% =0 

and 



(Q(f)J) 



s(v,v')M(v)F(x,v') 



t \ ( f(x,v') f{x,v) 



F(x,v') F(x,v) 



dv' dvdx 
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implying that f(x,v) = p(x)M(v) for some function p = p(x). Therefore Q(f) = and thus 
Qi(/) = 0. It is easily seen that for f(x,v) = p(x)M(v) the formula (|67|) reduces to 

Q\f) = -M(v)^ ■ [ (^p(x) - ^eP) ( | V | S1 _ \ v >\ S2 ) M (v>) dv>. 

We are ready if we prove that there is r > such that 

J(v) := / {\v\si - \v'\s 2 )M{v') dv' > 0, \v\ = r 



since in that case we deduce that 



u c (x) \p F (x) J u c (x) 

saying that p = Cp F for some constant C, or f(x,v) = Cp F {x)M{v) = CF(x,v). Recall that 
s\,S2 depend only on r = \v\,r' = \v'\ and 

si = — / o-(y/x) do, s 2 = — cosada, x = r 2 + ( r ') 2 ~ 2rr ' cos a - 

Therefore one gets 

J(r)= / (r- r r cos a)a(^/x)fj,(r') dadr' = - / / —a(^/x)p(r')dadr' 
Jr+ Jo 2 Jr+ Jo or 

m(r') 



where /u(r') = m/(2ir0) exp( ^r - )- Integrating with respect to r € M + yields 

/ J(r) dr = - / / M(r') / o^^x) dxdadr' > 0. 



Therefore there is r* > such that J(r*) > 0. ^ 

Another point to be discussed is related to the continuity equation. Since for any e > the 
solution f £ of the linear Boltzmann equation ([7|) satisfies a equation like 

d t [ f £ dv + div x [ vf £ dv = 



we should retrieve a similar one when integrating the second order approximation (|63p with 
respect to u £ t 2 . In the case of constant scattering cross section we obtain 

d t [ f dv + ediv x f (v A + v GD )f dv = e I Q l (f) dv = e-dw x [ —f dv 

JR2 JR 2 Jjg.2 T J R 2 UJ C 

and the continuity equation accounts for a new drift wqd = — ^ — ^ , orthogonal to the 
magnetic lines, due to collisions 

<9 t / / e dt; + ediv x / (v A + u G d + ^qd)/ £ dv = 0. 
Jr 2 jr 2 

In the general case, it is easily seen that 

I Q 1 (f £ )dv = -dw x f [ ( Sl -^L S2 ) M(v) — f £ (x,v')dv' dv = -div x f v QD f £ dv 
Jr 2 t Jm2j R 2\ \v I J U) c J R 2 

and therefore the collision drift is given by 

^QD = -- / (si - yjS 2 ) M(v') dv' — V —. (69) 

T Jr2\ \v\ J UJ c {x) 

The previous arguments provide the second order approximation model, as stated by Theorem 
12.31 For the sake of the presentation we assume here that the scattering cross section is 
constant. 
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Proof, (of Theorem 12. 3p Consider / the solution of the problem dtf = Q(f),t > with the 
initial condition /(0) = f m and g 1 the solution of the problem 

d t g 1 + (a-V x , v h 1 ) = Q(g 1 ), t>0, 5 1 (0)=0 

where h 1 = -B ■ V x , v f. Notice that Tg 1 = 0, (h 1 ) = and by the previous computations we 
know that 

d t f + a-V x , v f + Tf = Q(f), f 1 = g i + h 1 . 
We take also a smooth function f 2 such that 

d t f 1 + a-V x , v f 1 + Tf 2 = Q(f 1 ) 

which is possible thanks to the evolution equation satisfied by g . As before we write 

(d t + a- V x , v + e- 1 T)(/ + ef 1 + e 2 f 2 ) = Q(f + ef 1 + e 2 f 2 ) + e 2 (d t + a- W x , v - Q)(f 2 ). 

Combining to the equation satisfied by f £ we obtain, with the notation r £ = f £ — (f + ef 1 + 
e 2 f 2 ) 

(d t + a- V x , v + e- x TY = Q(r e ) + e 2 (Q-d t -a- V x , v )f 2 . 
Using the entropy inequality we deduce that, for some constant Ct 

\^V\^<Ct^V\\^, te[0,T], e > 

implying that 

\\r £ (t)\\ Ll <\\r £ (0)\\ Ll +TC T e 2 , te[0,T), e>0. (70) 

Notice that 

r £ (0) = f £ n -(f(0) +£ f 1 (0)+e 2 f 2 (0)) = f £ n -(r + eh 1 (0)) + O(e 2 ) 
= fin ~ (f m - eB ■ V x , v r) + 0(e 2 ) = 0(e 2 ) 

and therefore 

sup \\f(t)-(f(t)+ef\t))\\ L2 =0(e 2 ). (71) 

0<t<T F 

We also know (cf. (pTT j) . (p2~|) ) that / + ef 1 satisfies 

(d t - eC ■ V x , v )(f + ef 1 ) = (Q + eQ')(f + ef 1 ) - e 2 (C ■ V^/ 1 + Q 1 (f 1 )) 

which is a second order perturbation of the equation ([63]) (or dHJ)) satisfied by f £ , implying 
that 

(dt - eC ■ V X)V )r* = (Q + eQ v )(r £ ) + e 2 (C- V^f 1 + Q 1 (f 1 )), r £ = f - (f + ef 1 ). 
Using now the entropy inequality of Q + eQ 1 yields 

l_d 

2dt"' " l f ~ "■ 



\r £ f L 2 <C T e 2 \\r £ \\ a , te[0,T], e > 0. 



Taking into account that f £ (0) = f £ (0) - (f(0) + e^(0)) = / £ (0) - (/ in - eB ■ V x , v f in ) = 
we obtain 



F 



< TC T e 2 , t G [0,T], e > 



saying that 

sup ||f (t)-(f(t)+ef 1 (t))\\ L2 =0(e 2 ). (72) 

0<i<T b 

Our conclusion follows immediately from (|7ip . ()72p . | — | 
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A Proofs of Propositions 13.5 



3.6 



Proof, (of Proposition | 
1. We show that y — > (ipo(y), ...,ip m -i(y)) is a change of coordinates. Indeed, if y,y E M m 
verify ipi(y) = ipi(y), i E {0, 1, ...,m — 1}, then y,y belong to the same characteristic. Thus 
denoting by yo the discontinuity point of tpo on this characteristic, there are h,h E [0, T c (yo)), 
with h < h without loss of generality, such that y = Y(h;yo),y = Y(h;yo). Integrating 
(6° • V y ipo)(Y(s; yo)) = I(Y(s; yo)) = /(yo) between h and h we obtain 

o = MV) - My) = MY(h; yo)) - MY(h; yo)) = (h- h)i( yo ). 

Therefore h = h which implies y = y. We have shown that y E M m is uniquely determined 
by ipo(y), ...,ifj m -i(y). Indeed, y belongs to the characteristic associated to the invariants 
tpi(y), ...,ip m -i(y) and, if we denote by yo the discontinuity point of t/'o o n this characteristic, 
we have y = Y(s;yo), where the parameter s E [0,T c (y)) is determined by 

My) ~ Mvo) = f\b° • VyM(Y(r;yo)) dr = S I(y ). 



Finally, without loss of generality we suppose that *po(yo) = and thus ipo(y) £ [0, T c (yo)I(yo)) 
[0, [ipo] ) = [0,5"). Clearly the map y — > (V'o(y); V'm-i(y)) is a surjection between R m and 
[0, S) x D, which shows 1. 

2. Notice that V y ipo ^ spanjVy^i, V^Vm-i} since b° ■ V y ipo ^ and b° ■ V y Vi = ••• = 
6° • Vyip m -i = 0. Thus, for any i E {1, m — 1} there is a unique vector field b l such that 

/''•V.r, tfj, jE{0,l,...,m-l}, 

which proves 2. [— j 

Proof, (of Proposition 13. 6p 

Notice that for any y E M. m the function 

s -)• u(Y(s;y)) = w(ip ( Y ( s 'iy)),My)-^m-i(y)) 

is continuous on K. In particular this holds true for any discontinuity point yo of -00 • For any 

y = Y(s;y ), s E (0,T c (y )) we can write 

(6° • V y u)(y) = ^-{u(Y(s;y ))} = b° ■ V y ipo(y) ^o^(^o(y), ^lG/o)-, M-l(Vo)) 

= i(y) d tpo w (My)^-^m-i(y))- 

It remains to analyze the differentiability around the point yo- Without loss of generality we 
assume that / > 0. Taking s > one gets 

b° ■ V y u(yo) = lim -[w(i/j ( y (s; yo)), Myo), -,fpm-i(yo)) - w(Myo), -,M-i{yo))] 

s\0 S 

= lim -[w(i/] (yo) + si,Myo), -,M-i{yo)) - w(Myo), ■■■,M-i(yo))} 

s\0 S 

= i(yo)dip w + (Myo), -,M-i{yo)) (73) 
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where d^ Q w + stands for the right derivative of w with respect to i/jq . Taking now s < 0, using 
the S-periodicity of w with respect to i[>q, we obtain 

b° ■ V y u(y ) = lim -[w(ip (Y(T c + s;y )), ...,i[) m -i{yo)) - w(ip Q (y Q ), ^ m _i(y ))] 

s/-0 S 

= lim -[w(ip {yo) + (T c + s)I, ...,i> m -x(yo)) - w(ip (y ), ^ m -i(z/o))] 

s/OS 

= lim -[w(i>o(y ) + si, ...,ip m -i(y )) - w(ip (y ), ...,^ m -i(yo))] 

s/'O S 

= i(yo)d^ w- (^o(yo),-,-0m-i(yo)) (74) 

where d^ Q W- stands for the left derivative of w with respect to ipQ. Combining ([73]) . (|73|1 we 
deduce that w is differentiable with respect to tpo and at any point y £ M. m 

b° • V y u = I(y)d^ w(ip {y), ...,ip m -i(y)). 

Moreover, for any i £ {1, ...,m — 1} we have 

m— 1 

6 l • V y u = 9^w V ■ Vyipj = d^w(My), - .^m-id/)). y e M m . 

□ 
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